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ABSTRACT 

The actual value of the oxygen abundance of the metal-poor planetary nebula PNG135.9+55.9 has frequently been debated in the 
literature. We wanted to clarify the situation by making an improved abundance determination based on a study that includes both 
new accurate observations and new models. We made observations using the method of integral field spectroscopy with the PMAS 
instrument, and also used ultraviolet observations that were measured with HST-STIS. In our interpretation of the reduced and cali- 
brated spectrum we used for the first time, recent radiation hydrodynamic models, which were calculated with several setups of scaled 
values of mean Galactic disk planetary nebula metallicities. For evolved planetary nebulae, such as PNG135.9+55.9, it turns out that 
departures from thermal equilibrium can be significant, leading to much lower electron temperatures, hence weaker emission in colli- 
sionally excited lines. Based on our time-dependent hydrodynamic models and the observed emission line [Om] /15007, we found a 
very low oxygen content of about 1/80 of the mean Galactic disk value. This result is consistent with emission line measurements in 
the ultraviolet wavelength range. The C/O and Ne/O ratios are unusually high and similar to those of another halo object, BoBn-1. 

Key words. ISM: planetary nebulae: general - ISM: planetary nebulae: individual (BoBn-1, PN G135.9+55.9) - Hydrodynamics 



1. Introduction 

The stellar-like obj ect SBS 1150 +599 A from the Second 
Byurakan Survey dBalavanl 19971) has bee n spectroscopi- 
cally identified by Tovmassian et al.l (12001 1 hereafter T01; 
iGarnavich & Stanekll 19991) to be an old planetary nebula (PN) of 
the Galactic halo and was thereafter renamed PN G135.9+55.9. 
The same authors perform an abundance study based on pho- 
toionization models and come to the conclusion that this partic- 
ular object has the lowest oxygen abundance known so far for 
planetary nebulae, viz. below about 1/100 of the solar value. 

This object appears to be a challenge for any detailed spec- 
troscopic analysis since very few emission lines are detectable 
in the optical wavelength range. Because of a lack of suitable 
lines, it is impossible to make a direct plasma diagnostic, and 
it is also difficult to constrain photoionization models such that 
meaningful abundances will emerge. The analysis is, moreover, 
hampered since this object is faint, with a mean H/3 surface 
brightness of only ~ 5 x 10~ 16 ergcm~ 2 s -1 arcsec~ 2 (adopting a 
diameter of 5" and a total H/3 flux of 1.9xlO~ 14 ergcrrT 2 s _1 ; see 
Table [TJ, making it difficult to accurately measure lines close to 
the 1% level of H/3 (corresponding to ^2xlO~ 16 ergcrrT 2 s _1 ). 

These observational difficulties have led to a dispute about 
whether PN G135.9+55.9 is as extremely metal-deficit as T01 
claim. A critical point in this discussion is the determination 
of the stellar temperature from the ionization balance, using the 
line ratio of [Ne v] A3426 and [Ne in] A3869 in the nebula. The 
line strength of [Ne in] /13869 is uncertain, but decisive when 
fixing the effective temperature (T e ff) of the ionizing source. 
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iRicher et all d2002l hereafter R02) and lJacoby etall (120021 here- 
after J02) find a rather low r e ff 100 000 K, hence low oxygen 
abundances of less than 1/100 solar from the weak [O in] /15007 
line, the only oxygen line that is observed in the optical. 

iTovmassian et al.l ((2004, hereafter T04) conclude, by means 
of optical and ultraviolet (UV) spectra (FUSE), that the central 
ionizing source of PN G135.9+55.9 must be a very hot (T e g tss 
120 000 K) pre-white dwarf of rather low mass (« 0.55 M G ), 
which resides in a short-peri od binary system with a more mas- 
sive companion (P - 3 .92 h, [Narji wotzki et al]|2005l) . Recently 
ITovmassian et al.l 120071 based on X-ray observations, have es- 
timated that the temperature of this massive component is very 
high, r e ff=; 170 000 K, while the less massive optical-UV com- 
ponent is cooler, r e ff ^58 000 K. 

iPequignot & Tsamisl d2005l hereafter PT05), moreover, crit- 
ically examine the situation and come to the conclusion that all 
evidence favors a higher temperature of the central source, viz. 
r eff =* 130 000 K, and that the claimed strength of [Nem] ^3869 
most likely is wrong. The oxygen abundance would then be 
1/30-1/15 solar an d not as extreme as previously estimated. 
iJacobv et alj (120061 hereafter J06) also find a higher stellar tem- 
perature, using the much stronger UV emission lines of highly 
ionize d nitrogen and carbon to constrain T e g, IStasiriska et al.l 
d2005l hereafter S05) include UV lines in their study (likely us- 
ing the same HST-STIS data as J06) and conclude that the oxy- 
gen abundance is lower than what PT05 find, 1/130-1/40. 

Our aim was to better measure the nebular emission line 
spectrum, which is why we performed new observations using 
the integral-field spectrograph PMAS. With such a spectrum we 
could determine a reliable line strength, or an upper limit, of 
[Ne in] /13869. We begin in Sect. [2] with a description of our ob- 
servations and data processing, we then present our results in 
Sect. [3] We describe the physical setup of our radiation hydro- 
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dynamic models, how we estimate abundances and compare ob- 
servations with our models, in Sect. [4] In Sect. [5] we discuss non- 
equilibrium effects and compare our abundances with values of 
previous studies found in the literature. We close the paper with 
our conclusions in Sect. [6] 



2. Observations and data reduction 

Our observations were made with the 3.5 m telescope at Calar 
Alto using the lens ar ray (LARR) inte gral field unit (IFU) of 
the PMAS instrument dRoth et al.ll2005l) . The V600 grating was 
used to cover the spectral interval 3490-5150 A; at a dispersion 
of 0.81 ApixeL 1 and a resolving power of R = 1340. The wave- 
length range was chosen to cover the Balmer lines H/3 and Hy, 
as well as [Ora] AA 4959, 5007, Hen/14686, and [Nein] ,13869. 
The LARR IFU, furthermore, holds 16x16 separate fibers, where 
each fiber represents a spatial element on the sky. We used the 
0f.'5 sampling mode where every pointing with the IFU covers 
an area of 8" x 8" = 64" 2 on the sky. Any number of spatial 
elements can be co-added to create a final spectrum. 

Two 2700 s exposures, which were both centered on 
PN G 135. 9+55. 9, were taken at an airmass of 1.09-1. 12 on 2007 
February 13. Two 2700 s exposures were also taken at an airmass 
of 1.12-1.19 on 2007 February 14. These two latter exposures 
were offset by 5" E from the first two exposures. One additional 
1800 s exposure was taken at an airmass of 1.09 and an offset 
of 5" W in the second night. Weather conditions were less than 
optimal, and the seeing was 1'.'4-1'.'7. In addition to the science 
exposures continuum and arc lamp flatfields were taken, as well 
as spectrophotometric standard-star exposures of G191-B2B. In 
order to correct for a varying fiber-to-fiber transmission sky flats 
were taken at the beginning of the second night. 

Reducing the data we us ed the tool P3d_qnline that is a part 
of the PMAS P3d pipeline dBeckeri 120021: iRoth et al.ll2005b . At 
first the bias level was subtracted and cosmic-ray hits removed 
(cf. ISandin et al.l 120081) . Second, a trace mask was generated 
from an internal continuum calibration lamp exposure, identify- 
ing the location of each spectrum on the CCD along the direction 
of cross-dispersion. Third, a dispersion mask for wavelength cal- 
ibration was created using an arc-exposure. In order to minimize 
effects due to a significant flexure in the instrument all contin- 
uum and arc lamp exposures were taken within two hours of the 
respective science exposure. Fourth, a correction to fiber-to-fiber 
sensitivity variations was applied by dividing with an extracted 
and normalized sky flat-field exposure. In this process the data 
was changed from a CCD-based format to a row-stacked-spectra 
format. As a final step flux calibration was performed in iraf us- 
ing standard-star exposures. We did not correct for differential 
atmospheric refraction since we, due to poor spatial resolution 
and high seeing, have co-added the flux of 147 spatial elements. 
These elements fully cover the object with an area of36.75" 2 .In 
order to avoid the periodic line shift of stellar absorption lines of 
the central star (CS) we made a second spectrum where we did 
not add the nine spatial elements, which were the closest to the 
CS at ,Ih/? (we also used fewer elements on the outer boundary); 
this latter spectrum covers an area of 30. 75" 2 . 

H/3 is blended by the helium Pickering line He n ,14859, and 
the resulting line intensity is hereby about 5% too high (cf. J02; 
PT05). The other Balmer lines are likewise blended by other he- 
lium Pickering lines. Since we did not include the CS in our final 
spectrum we have not corrected for underlying stellar absorption 
as PT05 do. The discussion of the dereddening performed by 
R02 remains inconclusive, with the result of a somewhat uncer- 



tain and very low extinction towards PN G135.9+55.9. Also T04 
find a low interstellar extinction towards this object (E(B-V) 
0.04 ) that is mainly based on the Nhi column density they derive 
from FUSE spectra. Owing to the large differences of the line 
strengths found by different observers (see below), however, any 
(rather small) correction due to dereddening appears unimpor- 
tant at this point of our discussion. 

3. Results for PNG1 35.9+55.9 

We made all line fits using the IFU analysis p ackage ifsfit 
(Sandi n, in prep.) that was initially developed for lSandin et al.l 
d2008l) . We used a second order polynomial to fit the continuum 
and Gaussian curves to fit emission lines. Because of less favor- 
able weather conditions in the second night we only used the 
co-added spectrum of the two exposures of the first night, total- 
ing an exposure time of 5400 s, see Fig. Q] 

The intensities of all object lines we measured are presented 
in Table [U together with the values of R02, J02, and T01. We 
present both raw values, and values where we subtracted the con- 
tribution of Hen Pickering lines from the Balmer lines. Note 
that errors of our intensity measurements decrease for redder 
emission lines, reflecting the better signal-to-noise of the spec- 
trum in the redder region. Since the first night was not pho- 
tometric our value of the integrated flux of H/3 is very likely 
somewhat uncertain. Nevertheless, we calculated a spatially in- 
tegrated flux of H/3, using the spectrum where the CS is included, 
of 1.92x 10 _14 ergcirT 2 s _1 (Table Q]), and this value compares 
well with values given in the literature. R02 find values ranging 
from 1.0 .. . 2.6xlO~ 14 ergcirT 2 s" 1 , and J02 evaluate a total H/3 
flux of 1 .5 x 10~ 14 erg crrT 2 s _1 . While J02 correct their flux es- 
timate for the region that is not covered by the slit, the varying 
values of R02 seem to be a result of the width of the slit they use. 

We did not detect [Ne in] ,13869 (cf. the upper inset of 
Fig. [T). Since the signal-to-noise of our spectrum appears to be 
better than in previous studies, with five new measured lines in 
the nebula (see TableQ}, it is questionable if anyone has detected 
[Nein] ,13869. 

In order to calculate an upper detection limit of emission 
lines in our data we proceeded as follows. For a certain wave- 
length we assumed that we can detect a line with the maximum 
flux of cr above the continuum, where cr is the error of the mea- 
sured flux. We calculated the intensity / of this limiting line by 
integrating over a triangle of width 7 A, and / = 7/2 x cr. At 
the three wavelengths ,1 = 3790, 3860, and 5010 A we measured 

o-=6.17, 4.70, and 1. 46x1 0~ 17 erg cm _2 s _1 A~ . Dividing the re- 
sulting intensities with the intensity of H/3 we got the following 
limiting line ratios (xlOO): 1.1 (A = 3790 A), 0.85 (A = 3860 A), 
and 0.27 (A = 5010 A). Hence, we consider 0.01H/? as an upper 
limit of the line strength of [Ne in] ,13869. 



4. Data analysis using hydrodynamic models 

Consequences for estimates of abundances in the context of pho- 
toionization models and the assumption of thermal equilibrium 
are poorly known. In a study ba sed on hydrodynamical m odels 
using time-dependent io nization Schonbern er et al.l d2009l here- 
after Paper VII; also see lSchonberner et al.ll2005cl) demonstrate 
that metal-poor nebulae with low densities are prone to devia- 
tions from thermal equilibrium, because heating by photoioniza- 
tion is no more balanced by line cooling only. In extreme cases 
the electron temperature is also controlled by expansion cooling, 
which results in lower electron temperatures compared to the 
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Fig. 1. The co-added blue spectrum of PN G135.9+55.9 using 123 spatial elements, including the sky. The upper and lower panels 
show the wavelength ranges 3640-4350 A and 4320-5030 A, respectively. The insets in the upper and lower panels show a close-up 
of the spectrum where H8 and a potential [Ne ill] /13869, and [O in] AA4959, 5007 are found, respectively. Gray thick lines in the 
inset panels show the line fits. The wavelength ranges of the inset panels are indicated in the spectrum with horizontal lines. The 
positions of four emission lines of telluric origin are also indicated (Hg). The diffuse emission features at, e.g., A ^ 4510 A and 
A ^4980 A are also of telluric origin. For further details see Sect. [3] 



Table 1. Flux measurements of PN G135.9+55.9 



Emission line 


A [A] 


PMAS raw 


PMAS corr . 


R02 (SPM1) 


R02 (CFHT) 


J02 a 


T01 


/(case B) 


[Om] 


5006.84 


3.36 


(0.33) 


348 


(0.35) 


2.7 


(1.3) 


2.87 


(0.85) 


3 (1) 


3.1 




[Oiii] 


4958.92 


1.27 


(0.29) 


1.33 


(0.30) 
















w 


4861.32 


100.00 


(0.49) 


100.00 


(0.49) 


100.0 


(2.1) 


100.0 


(1.7) 


100 


100 


100.0 


He II 


4685.65 


78.72 


(0.62) 


82.11 


(0.64) 


76.1 


(2.3) 


78.6 


(1.5) 


77 (3) 


92 




Hen 


4541.59 


2.96 


(0.39) 


3.09 


(0.40) 
















Hy 


4340.45 


45.49 


(0.60) 


45.40 


(0.63) 


41.0 


(1.6) 


42.05 


(0.98) 


39 (3) 


56 


47.6 


He II 


4199.83 


2.48 


(0.59) 


2.59 


(0.61) 
















H5 


4101.74 


23.81 


(0.64) 


23.69 


(0.66) 


21.4 


(2.8) 


20.6 


(1.2) 


17 (2) 


30 


26.6 


H7 


3970.07 


9.70 


(0.73) 


10.11 


(0.76) 


10.0 


(2.5) 


6.11 


(0.64) 


4 (1) 


11 


16.4 


H8 


3889.06 


8.63 


(0.89) 


9.00 


(0.93) 






2.92 


(0.73) 






10.8 


[Nem] 


3868.80 


<1 




<1 








1.04 


(0.52) 


15 (7) 






H9 


3835.40 


8.22 


(0.99) 


8.6 


(1.0) 














7.5 


H10 


3797.91 


3.4 


(1.2) 


3.3 


(1.2) 














5.8 


F(H/3) b 




1.924 


(0.009) 






1.82 


(0.03) 


2.55 


(0.03) 


1.47 


1.19 





a The line ratios are corrected for an extinction value of E(B - V) = 0.03 and for a ~ 5% contribution of He n-lines to Balmer lines. 

b This is the total flux in emission measured for H/? in units of 10~ 14 erg cirT 2 s , systematic errors are not considered in the error estimate. 
Note. — Errors are given in parentheses. The emission line names and the rest wavelength /1 , and Case B line ratios for the Balmer lines to H/3, 
are given in Cols. 1, 2, and 9. We present our raw values in Col. 3 and values corrected for the contribution of He II Pickering lines to the Balmer 
lines in Col. 4. Columns 5-8 give the values presented by R02 (SPM1), R02 (CFHT), J02, and T01, respectively. 



(standard) equilibrium case by up to 30%, although ionization 
is still close to equilibrium. Diffe rences are insignifica nt when 
solar abundances are used instead dPerinotto et al.|[l998T) . 

In this section we combine our observed line strengths and 
UV-data obtained with HST-STIS (Jacoby, priv. comm.) with 
outcome of our time-dependent hydrodynamic models in order 
to estimate abundances. The UV-spectrum we used is publicly 
available and can be retrieved from the HST archive (see pro- 
posal 9466, PI: Garnavich); it is also published by J06 (see Fig. 1 
therein). 



At first we present our modeling approach in Sect. 14. II there- 
after our model sample in Sect. 14. 21 and a discussion of how we 
determine our abundances in Sect. 14.31 

4.1. Physical properties of our time-dependent models 

Our one-dimensional radiation hydrodyn amic (RHP) models of 
envel opes of PNe are described in detail in lPerinotto et al.l d!998l 
I2004I also see references therein). We emphasize that our models 
calculate ionization, recombination, heating, and cooling time- 
dependently at every time step. The cooling function is com- 



4 



C. Sandin et al.: The evolution of planetary nebulae VI. 



posed of the contribution of all considered ions, and for every 
individual ion up to 12 ionization stages are taken into account. 
Physical input parameters to the models include properties of 
the coupled CS model, element abundances, and the density and 
velocity structures of matter in the envelope. 

In calculating the models the wind of the CS was used as in- 
put at the inner boundary of the grid. The full model evolution 
was thereafter followed across the Hertzsprung-Russell diagram 
for about 15 000 years until (partly) recombination sets in close 
to the turn-around point (at an age of about lOOOOyr), and into 
subsequent stages of re-ionization due to advanced expansion. 
An important feature of our models is that we, at any time, can 
switch off all time-dependent terms and simultaneously fix the 
density structure and the radiation field. The models are there- 
after evolved until they settle into equilibrium, we are then able 
to study differences between dynamical and static models (also 
see Paper VII). We refer to these models as equilibrium models, 
in comparison to the dynamical models. We calculated surface 
brightnesses, emission line profiles and strengths of individual 
lines using a supple mentary code that is based on a version of 
lGesickietal.1dl996l) . 

4.2. Properties of our model sample 

A definite mass estimate for the ionizing source of 
PNG135.9+55.9 is still lacking. T04 identify the nucleus 
of PN G135.9+55.9 as a close binary, which hampers their 
analysis, using non-LTE model atmosphere spectra, to derive 
the mass of the ionizing star. They fit the photospheric Balmer 
lines of selected optical spectra to obtain the surface gravity 
\ogg. Assuming an effective temperature of T e ff = 120 000 K, 
and that the companion does not contribute to the optical 
emission, a comparison with stellar evolutionary tracks indicate 
a mass of 0.88 M that they infer as unrealistically high. 

T04 support the final estimate of 0.55-0.57 M Q considering 
only the Population II characteristics of the object (not being 
in the Galactic plane, having a high radial velocity and being 
metal-poor), and its relatively high kinematic age that the au- 
thors consider a rough estimate for the post-AGB age. In order 
to be an evolved CS, which still is in a pre-white dwarf stage, 
a kinematic age of fki n = 16 000 yr de manded a mass range that 
low (see Fig. 9 in T04). The claim of iTovmassian et aTl d2007l) 
that the binary core of PN G135.9+55.9 consists of one compo- 
nent of 0.565 M Q and T e s - 58 000 K, and a second component 
of 0.85 M Q and T e ff = 170 000 K, agrees poorly with their pre- 
vious results. For instance, the hydrogen and helium lines will 
hardly be fitted by an object where T e $ = 58 000 K. The ioniza- 
tion of the nebula is likewise inconsistent with an ionizing source 
where T efi ■ = 170 000 K, as PT05 show. 

PT05 set r e ff, the object distance and luminosity, and se- 
lect models using a range of masses (0. 583-0. 600M o ) based 
also on arguments using a kinematic age. Conclusions, which 
are based on kinematic ages (Rxxx/V yy y), are questionable when 
they are used as evolutionary ages. Errors emerge from uncer- 
tain distances, the confusion of matter (Doppler) velocities with 
structure (shock) velocities, inappropriate (i.e. unrelated) com- 
binations of R XXK a nd Vyw, and neglecting t he expansion his- 
tory of the object dSchonberner et all l2005ai also see Fig. 33 
in Paper VII). The subscripts (xxx and yyy) denote that var- 
ious combinations of R and V could be used in order to ob- 
tain a kinematic age. For R one could use the outermost radius 
(R ut), the rim radius, or the radius of any other substructure. For 
V one could, moreover, use corresponding differential changes 
(K xxx , although this quantity is rarely available), or alternatively 



Table 2. Mean Galactic disk element abundance distribution 

(Zgd) 



H 


He 


C 


N 


O 


Ne 


S 


CI 


Ar 


Z GD 12.00 


11.04 


8.89 


8.39 


8.65 


8.01 


7.04 


5.32 


6.46 



spectroscopic velocities such as the half-width half-maximum 
(HWHM), HW10%M of spatially unresolved profiles, or any ve- 
locity component of a decomposition of a spectrum of a spatially 
resolved profile. 

Because the properties and evolutionary history of 
PN G135.9+55.9 are uncertain we restricted our analysis 
to one stellar evolutionary track as input for our RHD mod- 
els, and used a post-AGB (CS) model of 0.595 M e . Using 
a similar approach as T04 we only considered one stellar 
component when evolving the nebula, and assumed that the 
influence from the companion is weak. The stellar wind from 
the C S is calculated according to iMarten & Schonbernerl 
d 199 lh . Its dependence on the metallicity (in this case C, N, 
and O) is approximately accounted for b y correction fac- 
tors; we used M oc Z 69 dVinket al.lfeOOll) and u ra oc Z° 13 
dLeitherer et all 1 19921) . whereby L wind = Q.5Mvi oc Z 095 . The 
CS radiates as a black body. For the AGB wind we, moreover, 
adopted a power-law density profile p oc r~" (a = 3-3.25, cf. 
Sect. 14.31 ) and a constant outflow velocity u = lOkms" 1 (cf. 



ISchonberner et al.1 l2005bl hereafter Paper II). We used the 
radial domain 4.0 X 10 14 < r < 2.8 x 10 18 cm. The models are 
normalized such that n = 10 5 cirT 3 at r — 3 x 10 16 cm. The 
abundances are based on scaled values of mean Galactic disk 
P Ne abundances (Zgd; th is abundance distribution is first quoted 
bv lPerinottoet~aTlll998l) for nine elements, see Table |2j except 
for carbon and nitrogen Zgd is close to solar. The abundance 
values e, are specified in (logarithmic) number fractions relative 
to hydrogen, i.e. e, = log (rii/nn) + 12. Abundance dis tributions, 
which are similar to Zgd, are used by for examp le, Perinottol 
d 19911). iKinpsburgh & Barlowl d 19941) . lExter et al.1 d2004l) . and 
iHvung etaL ( 20041) . We summarize all model parameters and 
properties in Table [4] Selecting metallicities we made use of 
the set of six models of Paper VII that we refined with four 
additional models (these additional models are marked with 
the prefix *): 3Z GD , Z GD , Z GD /3, Z GD /10, *Z GD /15, *Z GD /20, 
*Z GD /25, Zqd/30, *Z gd /60, andZ GD /100. 

In order to illustrate differences between dynamical and 
equilibrium models we present line strength ratios of three col- 
lisionally excited oxygen lines in the UV, optical, and infrared 
wavelength ranges in Fig. [2] Representing the model evolution 
the line ratios are shown as a function of the CS effective temper- 
ature T e g. The figure shows that lines of dynamical models are 
weaker than those of equilibrium models for T e g <: 50 000 K. 
The reason is that dynamical models have a lower electron tem- 
perature (cf. Sect. 3.1.3 and Fig. 16, both in Paper VII). We do 
not consider young PNe, where T e s < 50 000 K; in such objects 
the ionization equilibrium is disturbed for a short time, with the 
consequence that the ionization is somewhat overestimated in 
the equilibrium models. 

Discrepancies are larger with a lower metallicity, compare 
the thin lines with the thick lines in Fig. [2j and shorter wave- 
lengths, compare the dotted lines with the dashed lines. When 
we compare the line intensities at T e g ^ 130 000 K we see 
that, as expected, the difference is the largest for the UV- 
line O IV] /LI 1402+ 1405 that is up to 250% stronger in the 
equilibrium case of the Zgd/100 sequence. The optical line 
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Table 3. Observational vs. modeled line strengths 



- 1 1 1 1— 



Z QD /10 




150 



[1000K] 



Fig. 2. Comparison between emission line strengths of dy- 
namical (dyn) and equilibrium (eq) nebular models around a 
0.595 M central star evolving across the H-R diagram. The two 
shown model sequences use abundances Zgd/10 (thick lines) 
and Zqd/100 (thin lines). We show the dyn/eq line strength ra- 
tios (Ieq/Idyn) of three collisionally excited oxygen lines, viz. 
O IV] AAUQ2 + 1405 (dotted lines), [O ill] ,15007 (solid lines), 
and [O iv] A26 fim (dashed lines), as a function of T e ff. The evo- 
lution is only traced until maximum r e ff is reached in order to 
avoid confusion. See Fig. lA.3l for absolute values on the intensi- 
ties. For further details see Sect. 14.21 



[O in] A5QQ1 is stronger by about 70%. The infrared line is the 
least affected, it is up to about 30% stronger in thermal equilib- 
rium. 

Because differences in line strengths depend strongly on the 
input density and the velocity structure, these values should only 
be regarded as indicative. Nevertheless, non-equilibrium effects 
introduce an ambiguity to the determination of electron temper- 
atures, hence elemental abundances, which is neglected when 
using standard photoionization codes (cf. Sect.[5]i. 

4.3. Estimating abundances for PNG1 35.9+55.9 

A full abundance analysis of PN G135.9+55.9 is beyond the 
scope of this work, but as an example we determine abundance 
estimates that are based on our hydrodynamical models. Our 
model grid is not large enough to allow iteration of all param- 
eter dimensions. The primary goal of this study is instead to find 
a model that agrees reasonably with the observational quantities, 
and then use this best-match model to elaborate on the influence 
of non-equilibrium effects. 

We used the following four criteria to find such a best-match 
model: emission line strengths should match observed values, 
the model Ha emission line pr ofile should match t he observed 
line profile (such as Fig. 1 of iRicher et al.l 120031) . the model 
Ha surface-brightness distribution should match the observed 
distribution (such as Fig. 3 of R02), and both distance estimates 
from the H/3 flux of the object and from the corresponding appar- 
ent size, which is obtained from the surface-brightness distribu- 
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give the emission line names and the rest 
ues of PNG135.9+55.9 are given in Col. 
Corresponding values of dynamic (dyn) 
and their ratio, are given in Cols. 4-6 (cf. 



tion, should be in fair agreement. In addition to these properties 
that apply to the nebula, the visual magnitude of the central star 
is also indicative of the distance. However, since we have only 
used one single stellar track (with a stellar mass M = 0.595 Mq) 
such a visual magnitude is difficult to match properly. The emis- 
sion lines and their line strengths that we used in our study, are 
shown in Cols. 1-3 of Table[3] 

Since we have not been able to measure [Ne in] A 3869 we 
cannot determine T e ff using that line. Instead we proceeded as 
follows. Using the existing measurements of the two neon lines 
[Ne iv] AA 2422+2425 and [Ne v] A 3426 we determined a range 
of effective temperatures where these measurements match our 
model emission line strengths (see Fig. IA.4l i. The result is shown 
in Fig. [3^. The upper limit of the temperature range correspond- 
ing to the [Ne iv]-line is set by the maximum stellar temperature 
possible, in this case max(r e ff) = 146 870 K. Note that the extent 
of the two regions of the effective temperatures of [Ne iv] and 
[Nev] that do not overlap, depends on both observations and 
model initial conditions. We extrapolated both regions in order 
to find the nearest point of intersection (Z^ e , T e ff) that we then 
used as two of the initial values when iterating our models. From 
this procedure we found T e g = 138 000 K. 

For carbon, nitrogen and oxygen we evaluated line strengths 
at T e ff for all model sequences (Figs. lA~TT4A.3t : the variation of 
model line strengths with stellar effective temperature is mod- 
erate at r e ff = 138 000 K (see below). We show the results in 
Fig. |3j), together with power law fits in relevant regions. By a 
comparison with the observational data we then found a first set 
of object abundances. For sulfur, chlorine and argon we simply 
reduced the respective mean Galactic disk value using the mean 
scaling value that we found for carbon, nitrogen, oxygen, and 
neon. 

In a final step, during the creation of our first custom 
model, we adjusted the helium abundance. Figure IA.5I shows 
that model line strengths of helium are nearly independent of 
metal abundances. The same figure also shows that the varia- 
tion of each model sequence with the effective temperature is 
small for evolved objects (where T e g > 10 s K). We reduced the 
initial helium abundance by one third, as is suggested by the ra- 
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Table 4. Summary of our best-match model parameters 
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Fig. 3. For each model sequence panel a) shows the range 
of r e ff that corresponds to the observed line strengths of 
[Ne iv] AA2422, 2425 and [Ne V] ,13426, compare the plotted 
values with the neon line strengths in Fig. IA.4I The value of 
the effective temperature of our first custom model is indicated 
with an open circle (o; T s g = 138 000 K). In panel b) each sub- 
panel shows the line strength of all models (but 3Z) that are all 
evaluated at T e g = 138 000 K, for four different emission lines. 
The horizontal dashed lines indicate the observed values. In both 
panels error intervals are indicated with gray regions. All axes 
are logarithmic. For further details see Sect. 14. 31 



tio of model-to-observed line strength. With this set of initial 
abundances we then calculated additional custom models with 
modified abundances in order to converge on the observed line 
strengths. In this process we iterated element abundances, which 
were the closest to the observed values, accounting for the size 



Parameter 



PNG135.9+55.9 



Spectrum 
Stellar mass, M 
Model age, t 

Stellar effective temperature, r eff 
Stellar luminosity, L 
Central star wind, L^- mi = Q.SMv]^ 
AGB wind 



black body 
0.595 M 
8982 yr 
138 049 K 
2994 L 

M oc Z ° 69 , Va> oc Z 013 
p oc r~ 3 25 , v= 10kms~', 
n r =3xioi6cm=lxl0 5 cm- 3 



Abundances, e, =log (rc,/n H ) + 12: 



He 


10.88 


11.04 
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7.90 


8.89 
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7.47: 


8.39 
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6.74 


8.65 


Ne 


6.96 


8.01 


S 


[5.94] 


7.04 


CI 


[4.22] 


5.32 


Ar 


[5.36] 


6.46 


Distance, d 


18kpc 




Visual magnitude, m v 


19.5 mag 




Nebular density, (;i c ) 


65 cm" 3 




Nebular temperature, (r c > 


21 100 K 




Nebular H/8-luminosity, L(H/3) 


0.193 L G 




Model HWHM velocity, Vhwhm 


41.8kms-' 





Note. — The element abundances e, are used as input in the calculation 
of our radiation hydrodynamic models; the values of S, CI, and Ar are 
not fitted, but only scaled. Z G d denotes the mean abundance distribution 
in the Galactic disk (cf. Sect. l4.lt . 



of the error bars of the different lines. Setting the carbon abun- 
dance we first gave CIV/1/11548+ 1550 a stronger weight than 
C in] AA1906 + 1909, since its error is smaller (also see Fig.[3p). 

In the course of our iterative calculations we found that the 
model H/3 flux did not match the apparent size of the object at a 
common distance. In order to achieve a less extended model we 
modified the power law density distribution of the AGB wind by 
replacing a = 3, that is used in all other models of this study, 
with a = 3.25; keeping the density normalization as before (see 
Sect. 14. It . By this approach changes to resulting model emission 
line strengths should be small. In Paper II we show that a steeper 
density gradient results in faster expansion rates of the leading 
shock of the shell, i.e. the outer edge of the PN. Therefore our 
new choice of a may appear counterintuitive when forming a ge- 
ometrically smaller PN. The size of a PN, however, depends on 
the expansion velocity integrated over the entire expansion pe- 
riod. Early on a high circumstellar density at the inner edge of 
the model domain (due to our choice of a power law) prevents 
nebular matter from becoming rapidly ionized. The formation 
of a D-type ionization front is thereby delayed as the accelera- 
tion starts later with a steeper density gradient, and the overall 
PN lifetime thereby becomes shorter. In the case of a — 3.25 a 
higher expansion velocity cannot compensate for the simultane- 
ous shortening of the expansion period. 

The abundance distribution of our best-match model, and all 
relevant model properties, is given in Table|4j and resulting emis- 
sion line strengths are given in Col. 4 of Table [3] We also show 
model-to-observed line strength ratios in Fig. [4] This figure il- 
lustrates a weak dependence with effective temperature at values 
about T e s 138 000 K. Since most lines depend only weakly 
temperatures about this value, the precise value of T e g, for say 
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Fig. 4. This figure shows the evolution of emission line strengths of our best-match model, as a function of the effective temperature, 
T e ft. Using the ten observed emission lines of Table[3]we show model-to-observed line ratios (solid lines). The gray regions mark the 
corresponding error bars. Note that only upper limits were observed for O iv] AAIAQ2, 1405 and [Ne in] /I3869. For further details 
see Sect. 1431 



138 000 + 5000 K, is uncritical to the ionization structure. We 
could not achieve a simultaneous agreement for the two nitrogen 
lines. The high value of 426 for N v ^1238 + 1242 could not be 
reached with any of our models (Fig. IA.2l i. which is why the ni- 
trogen abundance of our best-match model should be considered 
approximate. Oiv] AA 1402+ 1405 can, furthermore, hardly be 
identified in the spectrum of J06. We consider the value of 37 a 
very conservative upper limit - compare with the value of 87 ±30 
for N iv] ,11486. A possible blending with Si iv ^11394+ 1403 
should not be excluded. 

In Fig. we show structural and kinematic properties of 
our best-match model; the CS has evolved to t — 8982 yr, L = 
2994 L , and r eff = 138 049 K (Table ©. The mean electron 
density is (n e ) = 65 crrT 3 . The density (Fig. |5^) shows a grad- 
ual decline with increasing radius. Neither the density nor the 
Ha surface-brightness structure (Fig. [5};) show a distinct double 
shell morphology. The same panel compares the observed data 
of R02 (see Fig. 3 therein) with the outcome of our model, us- 
ing a distance of d — 18.3 kpc. We show the real structure with a 
central dip that is due to the hot bubble, and the structure under 
the seeing conditions of the observations (l'/5). A comparison 
of our observed H/3 flux, F(H/3) = 1.924 x 10~ 14 ergcnT 2 s _1 , 
and the H/3-luminosity of the model, L(H/3) = 0.193 L Q , yields 
a distance of d = 17.95 kpc - that agrees well with the one of 



the surface-brightness structure. We adopt d — 18 kpc as a final 
value of the distance. In our models we use exclusively black 
bodies for the CSs. Applying the distance estimate of d « 18 kpc 
to the bolometric luminosity of the model, L = 2 994 L Q , yields 
my = 19.5 mag3 

Throughout the nebula the matter velocity gradient is posi- 
tive with increasing radius (Fig. |5^), reaching a maximum ve- 
locity of v - 65km s _1 . In the adjacent (radiative) shock layer 
at 12.5 SS r % 14.5xl0 17 cm the velocity is about constant. The 
simulated emission line pro file (Fig. |5}i) r esemb les the long slit 
Echelle spectral analysis of iRicher et al.1 (120031 see Fig. 1, slit 
7). Our one-dimensional model cannot be fully applied to this 
slightly non-spherical PN and its asymmetric line profiles. The 
observed HWHM velocity, Vhwhm = 42.5 kms~', is, however, 
well matched by our model with Vhwhm = 41.8 km s" 1 . 



1 R02 measure m v = 17.9 mag. Assuming this value instead our 
best-match model would shift to a corresponding distance of only 
d = 8.6 kpc. This discrepancy in our third distance estimate can be 
explained by our selected stellar mass (one single track of 0.595 M G ), 
which should also be iterated in order to achieve a better agreement. Of 
course, our model cannot reproduce the actual intensity if the nucleus 
really is a double degenerate. 
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Fig. 5. This figure shows the basic physical properties of our best fit model at the stellar parameters L = 2 985 L Q and T e g — 138 152 K, 
at an age of t — 8 988 yr. The four panels show: a) the radial structure of of the particle density (thick line) and the gas velocity (dotted 
line), b) radial structures of the electron temperature for the dynamical model (solid line) and the equilibrium model (dashed line), 
c) a comparison between the Hasurface-brightness distribution of the model at a assumed distance of 18.3 kpc and the observational 
data of R02 (filled/open squares = semi-major/semi-minor axis, see text for details) with the actual Hcr-image (dotted line) and 
an image that results with a seei ng of 1"5 (so l id line ), and d) a comparison between the Ha emission line profile of the model 
(solid line) and the observation of iRicher et al.l (I2003L see Fig. 1 and slit 7; open circles o). The simulated profile was additionally 
broadened by a Gaussian with FWHM= 26 km s _1 , in order to be consistent with their observation. The observed left wing is due to 
emission of He II /16560. For further details see Sect. 14.31 



In Fig. [5]5 we show the radial electron temperature structure. 
The temperature peak behind the outer shock does not contribute 
to the mean temperature due to low ion densities in that region. 
We will in the following subsection discuss the consequences 
of non-equilibrium conditions for the electron temperature and, 
consequently, for the strengths of collisionally excited lines. 

5. Discussion 

In order to study differences in the outcome of our time- 
dependent and static (equilibrium) models we should, ideally, 
also calculate a corresponding best-match equilibrium model 
by the same procedure we used to find the best-match dynam- 
ical model. As such an approach is extremely time-consuming 
we instead compare our outcome with literature values, which 
are all based on standard photoionization codes (these corre- 
spond to our equilibrium models). We present our abundances 
anew in Table [5] together with a compilation of literature val- 
ues, which are derived using observed emission lines. Errors 
of individual estimates are specified where such values are pro- 
vided. Differences between estimates and physical assumptions 
of different sources are large in general. Since we focus on un- 
derstanding general trends of values, and not on providing final 
abundances, we have not estimated errors of our values. This is 
also difficult to do with our models where abundances are input 
parameters, and not the outcome. 

Previous studies of PN G135.9+55.9 present improvements 
to different parts of the abundance analysis. PT05 provide a thor- 
ough model analysis, without making own observations, where 
they account for several physical issues, which were not ad- 
dressed previously. Notably, they study differences in models as- 
suming case B vs. non-case B photoionization, they use different 
sets of collisional recombination coefficients, and use a stellar 
atmosphere model of the CS, in addition to the commonly used 
black-body model. Due to lack of data they calibrate their mod- 



els using only observational data in the visual wavelength range, 
and therefore they cannot calculate precise values for the abun- 
dances of carbon and nitrogen. The main conclusion of PT05 is 
that the oxygen abundance of previous studies is too low. S05 
and J06, moreover, add UV lines, which are sampled with HST- 
STIS, to their linelist, and can thereby constrain the abundances 
of carbon and nitrogen better than PT05 (S05 also announce IR 
observations using SPITZER). S05 also argue that PT05 use too 
high abundances for carbon and nitrogen, and therefore have to 
use a higher oxygen abundance than is necessary; S05 find a 
lower value on the oxygen abundance than PT05, which is in 
better agreement with estimates of previous studies. 

In agreement with all previous studies, except for PT05, we 
found that the oxygen abundance of PN G135.9+55.9 is very 
low. It is difficult to make a meaningful, more detailed, compar- 
ison between our abundances and those of T01, J02, and R02 
since we use different stellar effective temperatures (mainly). 
PT05 are the latest authors who base their analysis on only op- 
tical emission lines. The very thorough analysis these authors 
make is found to be of small use as there is a strong disagreement 
between the predicted UV-line intensities of their best model 
(that agrees best with model Ml, cf. Table 3 in PT05) and ob- 
served UV line strengths (Table |3j; the difference is (with the 
exception of C in] /11908) in every case about a factor two. It is 
in this context worth mentioning that the electron temperature 
they adopt is about 9000 K higher than in our study, resulting 
in a different ionization structure (see below). J06 use a similar 
electron temperature as PT05 in their photoionization models, 
which is why their results should be affected to a similar degree. 
Furthermore, although S05 include UV-intensities in their anal- 
ysis they do not provide any temperatures at all and it is impos- 
sible to make a meaningful comparison with their abundances. 

Compared to the mean abundances of the Galactic disk our 
values for PNG135.9+55.9 are 1/1.45 (He), 1/9.8 (C), 1/8.3 (N), 
1/81 (O), and 1/1 1 (Ne). Compared to the total mean metallicity 
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Table 5. Literature compilation of abundances estimates for PN G135.9+55.9 
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11 Jacoby (priv. comm.) 

Note. — The table only includes estimates that are based on observed emission lines. Columns 1-4 specify the source reference, the wavelength 
range (O - optical, and U - UV), the stellar effective temperature (r cff ) and the mean electron temperature ((T e )) used in the study. Columns 5-9 
give element abundances using the same units as in Table[2] In Cols. 10-12 we also give abundance ratios relative to oxygen. Uncertainties are, 
were provided, given in parentheses. A colon indicates an uncertain value of our best-match model. The abundances of our best-match model 
are given in the row marked this work. In the last two rows we, for comp arison, give the mea n abundance distribution of the Galactic disk (Z GD ; 
Table|3 and the halo-PN BoBn-1 (the values of this object are taken from Howard et al. 1993). F° r further details see Sect.[U 
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Fig. 6. In panel a) we show the ratio between the emission line 
strength values of our best-match model (dyn) and the observed 
values (cf. Table [3] and Sect. |4.3l l. For all, but the two nitrogen 
lines, N V ^1238 + 1242 and N iv] /U486, the model values and 
observations agree within error bars. In panel b) we show emis- 
sion line strengths of the thermally relaxed best-match model 
(eq) relative to the dynamical model (dyn) for all line listed in 
Table [3] The vertical dashed lines in both panels indicate limits 
of different spectral domains. 



of the Galactic disk our value is 1/13. Our abundance estimates 
relative to oxygen are in better agreement with the values of an- 
other halo PN, BoBn-1. In this case our values of C/O, N/O, and 
Ne/O are 33-260% higher, although PNG135.9+55.9 is consid- 
erably more depleted of metals. 

In Fig. [6] we illustrate differences between line ratios of our 
best-match dynamic model and its corresponding equilibrium 
model, plotted as a function of wavelength; we present the same 



data in Table [3] Fig. shows that the observed line ratios are 
satisfactorily matched by the model (Sect. |4~3l l. Equilibrium-to- 
dynamic model line ratios are shown in Fig. [(J for all lines that 
we used with our best-match model (also including IR lines; at 
the assumed temperature T e s ^ 138 000 K, compare with Fig. [2] 
for the models Zgd/10 and Zgd/100). In this case differences 
can be larger than 100% in the EUV, about 50% in the UV, about 
20-30% in the optical wavelength range, and about 10-20% in 
the infrared wavelength range. A hint of the importance of accu- 
rate line ratios is given in Fig.[3p. If a model line ratio changes 
by a factor two it is easily seen that different abundances are re- 
quired to match the change. Although the non-linear response to 
the nebula of the full model is complex, which is why plots such 
as Fig.[3]3 are unsuitable when making quantitative estimates of 
abundances. 

Differences in line ratios between dynamical and equilib- 
rium models occur as a consequence of a different sensitivity 
of collisionally excited lines to the electron temperature. The 
mean temperature in the nebular region of the two models are 
(T e )dyn =21 100 K and (r e ) eq = 25 100 K, compare the two radial 
structures in Fig. \5)p, the difference is significant. The electron 
temperature of our evolved metal-poor models is determined by 
line cooling and expansion cooling. It is worth noting that al- 
though the oxygen abundance lies closer to Zgd/100 the mean 
model abundance is closer to Zgd/10, and it is this higher abun- 
dance that determines the physical structure of the object (see 
Fig.[2]and Paper VII, Figs. 15 and 16). In an observational study 
using the full wavelength range (EUV-infrared), where measure- 
ment errors are sufficiently small, there should be significant 
problems determining abundances of metal-poor objects using 
models that are unable to account for dynamical effects. 

6. Conclusions 

PN G135.9+55.9 is an extraordinary object as it is a metal-poor 
PN with the lowest oxygen abundance known. In order to clarify 
contradictory abundance determinations of PN G135.9+55.9 in 
the literature we made a new study of this object using a two-fold 
approach. At first we re-observed the nebula and could measure 
a more accurate spectrum in the visual wavelength range than 
has been done so far. Unlike previous observational studies we 
could only measure an upper limit of [Ne in] ,13869 of 0.01H/3, 



10 



C. Sandin et al.: The evolution of planetary nebulae VI. 



although we measured five new lines in the nebula. We therefore 
chose to base our estimate of the stellar effective temperature 
using supplementary UV-data of Jacoby (priv. comm.). 

In the second part of our study we used a newly calculated 
set of our radiation hydrodynamic models in order to determine 
abundances and study the influence of time-dependent effects. 
In this case such effects are found to be important, causing lower 
electron temperatures in the nebula. Resulting line strengths of 
dynamical models are lower than in the corresponding equilib- 
rium models (these models are relaxed after all time-dependent 
terms are set to zero). Consequently, different abundances are re- 
quired to match line strengths when using either approach. We 
found that it is only possible to make a self-consistent abun- 
dance determination using a dynamical model that is constrained 
using measurements in the entire wavelength range. Our final 
set of abundances of the five most abundant elements is: 1/1.45 
(He), 1/9.8 (C), 1/8.3 (N), 1/81 (O), and 1/11 (Ne), all with 
respect to the mean abundance of objects in the Galactic disk 
(Zgd)- The total metallicity is Zgd/13. Additionally, using our 
single 0.595 M evolutionary track we found an effective tem- 
perature of PN G135.9+55.9 of T eff 138 000 K, and a distance 
of d — 18kpc. This distance is about the double value assumed 
by PT05, but agrees well with the value of T04. The mean elec- 
tron temperature of our dynamical model is (r e ) = 21 000 K; this 
is 4000 K lower than the value of the corresponding equilibrium 
model. 

Although we believe that our approach provides a most sig- 
nificant improvement when determining abundances of metal- 
poor objects our modeling can be improved to provide more ac- 
curate values. At first one could consider to iterate more dimen- 
sions of the parameter space, such as e.g. the mass of the central 
star and properties of the AGB wind. Three additional sugges- 
tions for such improvements that are all considered by PT05 
are: using non-CaseB radiative transfer, replacing the black- 
body model of the central star with a model atmosphere, and us- 
ing improved collision rates. Observationally an accurate multi- 
wavelength study including the infrared wavelength range, such 
as is announced by S05, will help to constrain the models fur- 
ther. Last, but not least, it is important to clarify the parame- 
ters and evolutionary history of the ionizing central star(s) un- 
ambiguously. 
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Appendix A: Intensity evolution of our RHD models 

For each model sequence we mention in Sect. 14.21 we show the 
intensity evolution of all emission lines of Table|3]in Figs. [ATT} 
IA.5I Solid lines show dynamic models and dotted lines equi- 
librium models (for those sequences where they were calcu- 
lated). The abscissa is in every case the stellar effective tem- 
perature r e jf. Horizontal dashed lines mark observed values for 
PN G135.9+55.9, and gray shaded regions mark corresponding 
error intervals. Additionally we show the evolution of the two 
infrared lines, [Oiv] A26fim and [Ne v] Al4fim, despite a lack 
of currently existing observational data. 
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